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Abstract. The canonical front form Hamiltonian for non-Abelian SU(N) gauge theory in 3+1 
dimensions and in the light-cone gauge is mapped non-perturbatively on an effective Hamiltonian 
which acts only in the Fock space of a quark and an antiquark. Emphasis is put on the many-body 
aspects of gauge field theory, and it is shown explicitly how the higher Fock-space amplitudes can 
be retrieved self- consistently from solutions in the qq-space. The approach is based on the novel 
method of iterated resolvents and on discretized light-cone quantization driven to the continuum 
limit. It is free of the usual perturbative Tamm-Dancoff truncations in particle number and coupling 
constant and respects all symmetries of the Lagrangian including covariance and gauge invariance. 
Approximations are done to the non-truncated formalism. Together with vertex as opposed to Fock- 
space regularization, the method allows to apply the renormalization programme non-perturbatively 
to a Hamiltonian. The conventional QCD scale is found arising from regulating the transversal 
momenta. It conspires with additional mass scales to produce possibly confinement. 
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nonperturbative calculations 



^ 1 Introduction 



Over the past twenty years two fundamentally different 
pictures of hadrons have developed. One, the constituent 
quark model is closely related to experimental observation 
and phenomenology, see for example HJ^MQ and the ref- 
erences cited there. The other, quantum chromodynamics 
(QCD) is based on a covariant non-abelian quantum field 
theory. How can one reconcile such models with the need 
to understand hadronic structure from a covariant theory? 

Since "the Hamiltonian is of central importance in 
usual quantum mechanics" J5|, an intuitive approach for 
solving relativistic bound-state problems would be to solve 
the gauge-fixed Hamiltonian eigenvalue problem. But the 
presence of the square root operator in the equal-time 
Hamiltonian approach presents severe mathematical diffi- 
culties. Even if these problems could be solved, the eigen- 
solution is only determined in its rest system. Boosting 
such a wavefunction from the hadron's rest frame to a 
moving frame is as complex a problem as solving the 
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bound state problem itself. This reflects the contempo- 
rary conviction that the concept of a Hamiltonian is old- 
fashioned and littered with all kinds of almost intractable 
difficulties. Alternative procedures like those of Schwinger 
and Dyson or of Bethe and Salpeter are not easy to cope 
with either in practice, see for example ||. 

So, why does one not proceed with the only success- 
ful non-perturbative approach to gauge field theory, with 
lattice gauge theory and its recent quantitative pre- 
dictions §? 

It is the concept of a wavefunction which is so ap- 
pealing in an Hamiltonian approach, particularly if one 
chooses the front form of Hamiltonian dynamics Q . The 
light-cone wavefunctions encode the hadronic properties 
in terms of their quark and gluon degrees of freedom, and 
thus all hadronic properties can be derived from them. 
One can compute electro-magnetic and weak form factors 
rather directly from an overlap of light-cone wavefunctions 
pO[ , see also pH] , and the hadron and nuclear structure 
functions are comparatively simple probability distribu- 
tions. Lattice gauge theory has still a long way to go to 
formulate such dynamical aspects. 

As recently reviewed fLl), one has several reasons [|l2] 
why the front form of Hamiltonian dynamics is one of the 
very f ew candidates for getting wave functions, see also 
[OJlJ]. Particularly the simple vacuum and the simple 
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boost properties confront with the complicated vacuum 
and the complicated boosts in the conventional Hamil- 
tonian theory. Based on the recognition that "the only 
truly successful approach to bound states in field theory 
has been quantum electrodynamics (QED) with its com- 
bination of non-relativistic quantum mechanics to han- 
dle bound states and perturbation theory to handle rel- 
ativistic effects" Jl9], Wilson and collaborators Cq] have 
proposed a scheme in which one presumes a potential for 
the bound-states and handles the relativistic effects by 
structures imposed by the needs of renormalization. This 
scheme continues to develop |l(|. The only available nu- 
merical examples |l7],|l8|,[l9| however violate each and ev- 
ery symmetry of the Lagrangian. 

The aim and motivation of the present work is sim- 
ilar to One aims at finding an effective interaction 
between quarks in analogy to the Coulomb potential as a 
crude and zero-order approximation to a given Lagrangian 
which can be used as a starting point for more refined con- 
siderations within a Hamiltonian approach. 

But the problems one faces with a Hamiltonian are 
stupendous. One has to deal with many difficult aspects 
of the theory simultaneously: confinement, vacuum struc- 
ture, spontaneous breaking of chiral symmetry (for mass- 
less quarks), and describing a relativistic many-body sys- 
tem with unbounded particle number. The problem in 
non-Abelian gauge theory is compounded not only by the 
physics of confinement, but also by the fact that the wave 
function of a composite of relativistic constituents has to 
describe systems of an arbitrary number of quanta with 
arbitrary momenta and helicities. 

For example, the vacuum in the front form is not truly 
simple [ pOpl ] but still simpler than in the instant form. 
The problem can at least be formulated in terms of the 
'zero modes', the space-like constant field components de- 
fined in a finite spatial volume. The original conjecture 
that zero modes represent long range aspects and thus 
confinement p2| , however has not materialized p3| . Zero 
modes are important for a 'theory of the vacuum', but 
probably less for the massive part of the spectrum. In 
the present work zero modes will therefore be suppressed 
without any good argument other than simplicity. With 
similar arguments all aspects of chiral symmetry (break- 
ing) are disregarded. 

One addresses thus to find the bound-state structure of 
the light-cone Hamiltonian in the light-cone gauge A + = 
0. This gauge is kind of natural to the light cone, and 
it is physical since the gluons have only two transverse 
degrees of freedom. Simple considerations show why other 
than in 1+1 dimensions one should not address to solve 
the full Hamiltonian |Tl| in physical space-time. In the 
first place one should address to develop a well-defined 
effective interaction, such one as proposed for example by 
Tamm and independently by Dancoff |^| , in their 
paradigmatic treatment of the Yukawa model. But if one 
does so and adapts the Tamm-Dancoff approach to the 
light-cone [^6| , one finds no trace of a possible confinement 
in these 'mesons'. This is rather disturbing since the same 
approach applied to positronium yields the Bohr aspects 



of the spectrum including the correct fine and hyperfine 
splittings. 

The Tamm-Dancoff approach will be re-analyzed in 
Sect. H The problems are generic, since for the purpose of 
practical calculations one must truncate at the lowest non- 
trivial order of a perturbative series. If one tries to correct 
for that one gets into all kinds of difficulties including the 
question how one should resume perturbative series to all 
orders without double counting. The problem is a severe 
one since one is forced to mistreat precisely all those multi- 
particle configurations which one thinks are important for 
confinement in a field theory. Quite naturally, they come 
in certain higher powers of the coupling constant and are 
thrown away when one truncates the perturbative series. 
In a way that part of the problem is similar to throwing 
away many-body amplitudes in conventional many-body 
problems. The problem in field theory are accentuated by 
the fact that the terms of higher order in the coupling 
constants are to be multiplied with very large numbers. 

The shortcomings of the Tamm-Dancoff approach are 
overcome here by the method of iterated resolvents |27| , |28| J^ 
This novel method is inspired by the Tamm-Dancoff ap- 
proach and requires the development of a considerable for- 
mal apparatus. As to be shown, it allows for a systematic 
discussion of the many-body effects in a field theory with- 
out truncating the expressions at a finite order of pertur- 
bation theory, at any time. With only a few well localizable 
assumptions one ends up with integral equations, whose 
structure is similar what was solved already in the past 
[ p6|]30| , |3l| ] . Most importantly one arrives at the conclu- 
sion that the effective interaction has essentially only two 
contributions: a flavor-conserving and a flavor-changing 
part. In a diagrammatical representation both look like 
low order perturbative graphs. In an approximative sense 
they even are identical with perturbative graphs, with one 
important difference however: All genuine many-body ef- 
fects reside in the effective coupling constant and generate 
a dependence on the momentum transfer across the ver- 
tex. They generate new interactions and possibly generate 
confinement. 



2 The light-cone Hamiltonian (matrix) 

The advantages and challenges of quantizing field partic- 
ularly gauge field theory 'on the light cone' and of solving 
the bound state problem for quantum chromodynamics 
(QCD) have been reviewed recently pjj]. It should suffice 
therefore to recall here the most elementary facets and to 
shape some notational aspects in short. 

In canonical field theory the four commuting compo- 
nents of the energy-momentum four-vector P v are con- 
stants of the motion. In the front form of Hamiltonian 
dynami cs M they are denoted by P v = (P + , Pj_, P~), see 
also p"3| , [ll[ ]. Its three spatial components = (P^P 2 ) 
and P~^ do not depend on the interaction. The fourth com- 
ponent P~ = 2P + depends on the interaction and is a 
complicated non-diagonal operator. It propagates the sys- 
tem in light-cone time x + = x°+x 3 , i.e. i-^rW) = P+l^), 



H.C. Pauli: On confinement in a light-cone Hamiltonian for QCD 



3 



Table 1. The Fock space and the Hamiltonian matrix H' = T + V for a meson at fixed value of K — 4. — See discussion 
in the text. The diagonal blocks are denoted by T. Most of the block matrices are zero matrices, marked by a dot (■). The block 
matrices marked by V are potentiall non-zero due to the vertex interaction. 
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and is the front-form Hamiltonian proper || . The contrac- 
tion of P M is the operator of invariant mass-squared, 

P^P„ = P+P- - Pf =H LC = H . (1) 

It is a Lorentz scalar and referred to somewhat improp- 
erly, but conveniently, as the 'light-cone Hamiltonian -Hlc' 
|l3| , or shortly H. Its matrix elements (and eigenvalues) 
somewhat unusually carry the dimension of an invariant- 
mass-squared. 

In this work one aims at a representation in which H 
is diagonal, 

H\#i) =E i \W i ). (2) 

One wants to calculate the bound-state spectrum Ei to- 
gether with the corresponding eigenfunctions directly 
from the gauge field QCD-Lagrangian. 

A convenient Hilbert space for the energy-momentum 
operators P v is the Fock space. The Fock space is the 
complete set of all possible Fock states 

|* B >= blX--- «---|0). (3) 

The creation and destruction operators obey the famil- 
iar (anti-)commutation relations. Each particle has four- 
momentum p^ = (p + ,p±,p~) and sits on its mass-shell 
p^p^ — m 2 with (light-cone) energy p~ — (m 2 + p 2 )/p + . 
Each quark state is specified by the six quantum numbers 
q = (p + , pi, A; c, /), i.e. the three spatial momenta, helic- 
ity, color and flavor, respectively. Correspondingly, gluons 
are specified by five quantum numbers q — (p + ,px, A; a) 
with the glue index a. The three spatial components 
and P + are diagonal operators in Fock space representa- 
tion with eigenvalues 

P+ = ]Tp+ , and Pj_ = £>J-)n . (4) 

n n 

The matrix elements of P~ (or of H) in Fock-space rep- 
resentation are given elsewhere fll|| . 



The eigenvalue equation in Eq.(g) stands for an infi- 
nite set of coupled integral equations which are extremely 
difficult to handle, see It is useful to convert it to the 
much more transparent case of a finite set of coupled ma- 
trix equations, namely by the technical trick of imposing 
periodic boundary conditions (DLCQ, |$^]). Introducing 
a box size L as a finite and additional length parameter, 
however, can be at most an intermediate step. Latest at 
the end of the calculations, it must be removed by a lim- 
iting procedure like L — ► 00, K — ► 00, but K/L finite, 
since only the continuum can be the full, covariant theory. 

Why is this set finite? — The longitudinal light-cone 
momentum p + is a positive number. For periodic bound- 
ary conditions the lowest possible value is (p + ) m i n = tv/L. 
(Zero modes with p + = arc disregarded here, as men- 
tioned.) Consequently, any total momentum P + = Kir/L 
can be distributed over at most K bosons, or over K 
fcrmion pairs since these are subjected to anti-periodic 
boundary conditions. As illustrated in Table [I] for the Fock 
space of a meson, the harmonic resolution K fl32| governs 
the number of Fock space sectors. 

The lowest possible value K = 1 allows only for one 
Fock-space sector with a single gg-pair (a single gluon can- 
not be in a color singlet). For K = 2, the Fock space 
contains in addition two gluons, a gg-pair plus a gluon, 
and two gg-pairs. For K = 4 the Fock space contains at 
most 8 particles. One can label the Fock space sectors 
according to the the quark-gluon content, or one can enu- 
merate them, which is less transparent but more simple. 
In Table |l| the Fock-space sectors for K < 4 are enu- 
merated n = 1, . . . , 13. With increasing K more Fock- 
space sectors are added. Their total number grows like 
Njc = (K + 1)(K + 2)/2 — 2. 

Once the Fock space is defined, one can calculate the 
Hamiltonian matrix. In jll] H is divided into three struc- 
turally different parts: 

H = T + V + W. (5) 
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The kinetic energy T is independent of the coupling con- 
stant g. It is diagonal in Fock-space representation, and 
contributes only to the diagonal blocks in Table [j]. These 
diagonal blocks carry only diagonal matrix elements which 
are given by Eq. (||) . — The vertex interaction V is the rel- 
ativistic interaction per se and linear in g. The odd number 
of creation and destruction operators prevents diagonal 
matrix elements. Potentially, the vertex interaction has 
non-zero matrix elements only between Fock-space sectors 
whose particle numbers differ by 3 or by 1 . Those differing 
by 3 are the typical vacuum fluctuation diagrams in which 
a fermion pair is created simultaneously with a gluon. On 
the light-cone, they vanish kinematically because of lon- 
gitudinal momentum conservation: The vacuum does not 
fluctuate, see for example |n]]. In consequence, the ver- 
tex interaction on the light cone can change the particle 
number only by 1. The corresponding blocks are marked 
by a V in Table [l| — The instantaneous interaction W 
arises as a consequence of working in the light-cone gauge. 
Potentially, it has non-zero matrix elements between Fock 
space sectors whose particle number differs by or 2. In 
the following discussion W will be left out because of a 
much more transparent formalism. The impact of W will 
be restored by a simple trick towards the end of Section [|. 

Table |l| highlights some of the main issues of a Hamil- 
tonian approach. It demonstrates that the Hamiltonian 
matrix is very sparse: most of the block matrices between 
the sectors are plain zero matrices. Very much like in the 
non-relativistic case with its typical pair-interactions, the 
Hamiltonian cannot connect all Fock-space sectors due to 
the selection rules imposed by the creation and destruction 
operators. Depending on the arrangement of the sectors, 
the Hamiltonian matrix has a (tri-diagonal) band struc- 
ture. The table illustrates further that the Hamiltonian 
'on the light cone' is separable into a kinetic part and 
the interaction, again like in conventional non-relativistic 
many-body problems. Therefore it should be approach- 
able like that, namely by introducing some complete and 
denumerable set of functions |n), in terms of which one 
can calculate the 'Hamiltonian matrix' (n\H\m). Its diag- 
onalization, 

Y,(*jW)(n'\H\n)(n\# i ) = E i 5 ij , (6) 

is equivalent to solving Eq. (0) . The columns of the unitary 
transformation matrix (n|^) are the 'Fock-space ampli- 
tudes' or 'wave functions'. The eigenstates are complicated 
superpositions of them, i.e. = J2 n \ n ) ( n l^)- 

But the analogue with non-relativistic Hamiltonian is 
superficial. There, the diagonal blocks (n\H\n) have off- 
diagonal matrix elements (potential energies) and one can 
approach the problem by successive truncation. In (gauge) 
field theory, the diagonal blocks contain only the (diago- 
nal) kinetic energies. Possible binding effects most come 
from the virtual scattering to the higher sectors. Trun- 
cating the matrix as done so often in many-body the- 
ory prevents such virtual scattering and potentially vi- 
olates gauge invariance. Moreover, since the theory ex- 
poses divergies which can be controlled by a cut-off or 



regulator scale A, all eigenvalues will depend on it: E- L = 
Ei(g, mf] A). Since this is unphysical, one must require for 
all of them that they are independent of A, i.e. 

± El (g,m f ;A) = 0. (7) 

Non-perturbative renormalization has not yet been ap- 
plied to a Hamiltonian, in practice. 

The bottle neck of any Hamiltonian approach as illus- 
trated in Table |f| as well: The dimension of the Hamil- 
tonian matrix increases exponentially. To give an exam- 
ple, suppose the regularization procedure allows for 10 
discrete momentum states in each direction, i.e. in the 
one longitudinal and the two transversal directions of kj_ . 
A particle has then roughly 10 3 degrees of freedom. A 
Fock-space sector with n constituent particles has thus 
10 n_1 different Fock states, since they are subject to the 
constraints, Eq.([|). Sector 13 alone with its 8 particles 
has thus the dimension of about 10 21 Fock states. The 
Hamiltonian method applied to gauge theory therefore 
faces a formidable matrix diagonalization problem. Sooner 
or later, the matrix dimension exceeds imagination, and 
other than in 1+1 dimensions one has to develop new tools 
by matter of principle. 

Aiming at an effective interaction between a quark 
and an antiquark, different novel methods have been pro- 
posed recently. Glazek and Wilson |Q propose to pre- 
diagonalize the Hamiltonian approximatively but analyt- 
ically such that its band width becomes narrower. The 
work on an iterative procedure is still going on [ [l6| . First 
applications to heavy mesons [fl8| have been carried out, 
but it is unclear how one can correct for the admitted vi- 
olation of every possible symmetry. Wegner |33|] has pro- 
posed an analytic unitary transform which leads to Hamil- 
tonian flow equations. Applications to realistic models [|4| 
are promising. Preliminary studies on QED [p5[ are avail- 
able. In parallel and partially prior to these works, the 
method of iterated resolvents p7[[2^ , |29| has been proposed. 
Its consequences are worked out in the sequel. 



3 Fock-space and vertex regularization 

Before proceeding with the theory of effective interactions, 
the regularization of the theory is specified next in order 
to face a well defined and finite Hamiltonian matrix. 

As mentioned, the finite number of Fock-space sectors 
is a consequence of the positivity of the longitudinal light- 
cone momentum p + . The transversal momenta pj_ can 
take either sign, and the number of Fock states within 
each sector can be arbitrarily large. In order to face a 
finite dimensional Hamiltonian matrix one must have a 
finite number of Fock states, and this is achieved by Fock 
space regularization: Following Lepage and Brodsky flipf , 
a Fock state with n particles is included only if its kinetic 
energy T, 

n / 2 I l , 2 \ 

T n ^(p 1 +p 2 + ... Pn f = Y,( !IL± ^ L ) • (8) 
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does not exceed a certain threshold. T can be interpreted 
as the free invariant mass (squared) of the Fock state. 
The lowest possible value of T n is taken when all par- 
ticles are at rest relative to each other, i.e. (T n ) min = 

(mi + m 2 + . . . m n ) . This frozen invariant mass should 
be removed from the cut-off 



Note that the semi-circle intersects the z-axis at z\ = q 
and z 2 — 1 — € u , with 



b-c 



Q = 



and 



(14) 



E 



(mi + m 2 + . . . m n ) < A 



FS 



(9) 



Since x and kj_ are the usual momentum fractions and 
intrinsic transversal momenta, respectively, the regular- 
ization is frame-independent . The mass scale /Ifs is a 
Lorentz scalar and one of the parameters of the theory. 

However, it was not realized in the past |i~3] ], that 
Fock-space regularization is almost irrelevant in the con- 
tinuum theory. Vertex regularization is a better alterna- 
tive. - At each vertex, a particle with four-momentum 
pV" is scattered into two particles with respective four- 
momentum Pi and P2- Parameterizing the momenta as 
p^ = (p+,p ± ,p~), p 1 ^ = (zp + ,zp± + l±,Pi) and p£ = 
((1 — z)p + , (1 — z) p± — 1l,P2~), tne (vertex) matrix ele- 
ment |p| is proportional to l_i/z. It tends to diverge for 
l± — ► oo and/or z — > 0. In order to avoid potential sin- 
gularities one can regulate the interaction by setting the 
matrix element to zero if the off-shell mass (pi + P2) 2 ex- 



For sufficiently large A and equal masses (mi = m 2 = m), 
the limits become ei = e u = (m/A) 2 . 

The scale parameter A regulates potential transversal 
divergences (l± — > 00), while potential longitudinal singu- 
larities (z — > 0) are regulated by the single particle mass 
m. If all particles are endorsed with a small additional 
regulator mass according to 



reg 



(15) 



ceeds a certain scale A. The condition 

ipi +P2) 2 - (mi +m 2 ) 2 < A 2 



(10) 



will be referred to as vertex regularization. The off-shell 
mass 

(Pi + P2) 2 - (mi + m 2 ) 2 = 



one ensures that the point z — is never inside the semi- 
circle given by Eq. (|l3|) , even not for the massless gluons, or 
for the quarks in the limit rrif —> 0. Vertex regularization 
regulates then all divergences on the light cone. 



4 The method of iterated resolvents 

Effective interactions are a well known tool in many-body 
physics [ j36| . In field theory the method is known as the 
Tamm-Dancoff-approach, applied first by Tamm and 
by Dancoff pEp . Let us review it in short. 

As explained above, the Hamiltonian matrix can be 
understood as a matrix of block matrices, whose rows and 
columns are enumerated by % = 1, 2, . . .N in accord with 
the Fock-space sectors in Tabic [l]. Correspondingly, one 
can rewrite Eq.(0) as a block matrix equation: 



z(l-z) 



N 

(V + (mi + m 2 ) 2 (z z) 2 ) , (11) (*l*> = E 



for all i = 1,2, 



governs how much the particles can go off their equilibrium 
values l± = and z = mi/(mi + m 2 ). For A — > 0, the 
phase space is reduced to a point, and consequently the 
interaction is reduced to zero: The Hamiltonian matrix 
(or the integral equation) is diagonal. Irrespective of the 
matrix dimension governed by A-ps, the spectrum of the 
interacting theory is identical with the free theory. 

Vertex regularization Eq. ([!(]) is realized by the cut-off 
function 0, i.e. 



1, for <l'l <l' 2 A {z), ei<z<l-e v 
0, otherwise. 



The limiting momentum l\(z) describes a semi-circle in 
the appropriate units, 



l\{z) = (A 2 + (mi + m 2 ) 2 ) 
yl 2 + 2mi(mi + m 2 ) 



The rows and columns of the matrix can always be split 
into two parts. One speaks of the P-space P — Y^j=i li) 01 
with 1 < n < N, and of the rest, the Q-space Q = 1 — P. 
Eq. (^6|) can then be rewritten conveniently as a 2 x 2 block 
matrix equation 

(P\H\P) (P\W) + (P\H\Q) (QW) = E (P\W), (17) 
(Q\H\P) (P\9) + (Q\H\Q) (Q\9) - E (Q\V). (18) 

One rewrites the second equation as (Q \ E—H \ Q) {QYP) = 
(Q\H\P) (P\'l r ), and observes that the quadratic matrix 
(Q\E — H\Q) could be inverted to express the Q-space in 
terms of the P-space wavefunction. But here is a problem: 
The eigenvalue E is unknown at this point. One therefore 
solves first an other problem: One introduces the starting 
point energy lj as a redundant parameter at disposal, and 
defines the Q-space resolvent as the inverse of the block 
matrix (Q\ui — H\Q). The Q-space function becomes then 



ith b = 



and 



A 2 + (mi + m 2 ) 2 

A 2 A 2 



4mim 2 



A 2 + (mi + m 2 ) 2 A 2 + (mi + m 2 ) 2 



(13) 



(Q\W(w)) = G q (lo)(Q\H\P)(P\*}, 
GqH = {Q\u-E\QY 



(19) 
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If there is no danger of confusion, the argument in the 
resolvents G(u>) will often be dropped in the sequel. In- 
serting it into Eq.(0) defines the effective Hamiltonian 

H eS (uj) = H + H\Q) G q {uj) (Q\H . (20) 

By construction it acts only in the P-space 

ifeffM|P) (P|#*(w)) = E k (w) |!Pit(w)) . (21) 

In addition to the original Hamiltonian in the P-space, the 
effective Hamiltonian acquires a piece where the system is 
scattered virtually into the higher sectors represented by 
the Q-space, propagating there (Gq) by impact of the true 
interaction, and scattered back into the P-space. Every 
value of lu defines a different Hamiltonian and a different 
spectrum. Varying lj one generates a set of energy func- 
tions Ek{uj). Whenever one finds a solution to the fix-point 
equation 

E k {u)=u), (22) 

one has found one of the true eigenvalues and eigenfunc- 
tions of H, by construction. 

One should emphasize that one can find all eigenval- 
ues of of the full Hamiltonian H, irrespective of how small 
one chooses the P-space. Explicit examples for that can 
be found in p9[ . It looks as if one has mapped a diffi- 
cult problem, the diagonalization of a huge matrix onto 
a simpler problem, the diagonalization of a much smaller 
matrix. But this true only in a restricted sense. One has to 
invert a matrix. The numerical inversion of a matrix takes 
about the same effort as its diagonalization. In addition, 
one has to vary ui and solve the fix-point equation (^2|). 
The numerical work is thus rather larger than smaller as 
compared to a direct diagonalization. 

The advantage of working with an effective interaction 
is of analytical nature to the extent that resolvents can be 
approximated systematically. The two resolvents 



GoM 



1 



(Q\lu-T -U\Q)' 
1 

{Q\u>-T\QY 



and 



(23) 



defined once with and once without the non-diagonal in- 
teraction U, are identically related by Gq = Go + GqUGq, 
or by the infinite series of perturbation theory. 



Gq = Gq + GqU Gq + GqUGqUGq + 



(24) 



The point is that the kinetic energy T is a diagonal opera- 
tor which can be trivially inverted to get the unperturbed 
resolvent G . 

In practice, Tamm and Dancoff |24|,^5| have restricted 
themselves to the first non-trivial order, and also the re- 
cent applications to the light-cone Hamiltonian have not 
gone beyond that j2^|nj]. This is unsatisfactory, since it 
destroys Lorentz and gauge invariance. Even worse, if one 
identifies uj with the eigenvalue (as one should), the ef- 
fective interaction develops a non-integrable singularity 
|m| , [25| . In the front form work (2^^] it was argued why 



the so called uj = uj* trick removes this singularity and 
approximatively restores gauge invariance. In the essence, 
one replaces the eigenvalue ui by the average kinetic energy 
in the P-space. 

The Tamm-Dancoff approach can be interpreted as the 
reduction of a block matrix dimension from 2 to 1, particu- 
larly the step from Eqs.([l7|jl|) to Eq.(|2l]). But there is no 
need for identifying the P-space with the lowest sector. In 
the sequel one chooses the last sector as the Q-space: The 
same steps as above reduce then the block matrix dimen- 
sion from N to N — 1. The effective interaction acts in the 
now smaller space. This procedure can be iterated. The 
disadvantage is that one deals with 'resolvents of resol- 
vents', or with iterated resolvents. The advantage is that 
the zero block matrices simplify the algorithm consider- 
ably. In the Tamm-Dancoff approach they cannot be made 
use of. Ultimately, one arrives at block matrix dimension 1 
where the procedure stops: The effective interaction in the 
Fock-space sector with only one quark and one antiquark 
is defined unambiguously [ p9| . 

Suppose, in the course of this reduction, one has ar- 
rived at block matrix dimension n, with 1 < n < N. De- 
note the corresponding effective interaction H n {uj). Since 
one starts from the full Hamiltonian in the last sector N, 
one has to convene that ipv = H. The eigenvalue problem 
corresponding to Eq.(^l|) reads then 

n 

5>|if»M|j>0W")) = E{d) (iMu,)), (25) 

3=1 

for i = 1,2, ... ,n. Observe that i and j refer here to sector 
numbers. Now, in analogy to Eq.([L9h, define 



(n|tf(w))= G„M^(n|P„M|j) (j|#H>, 

(n\u) - H„(u>)\n) 
The effective interaction in the (n— l)-space becomes then 



H, 



r „_i(u;) = H n (u) + H n (w)G n (u)H n (u) . (27) 

This holds for every block matrix element (i\H n -i(cj)\j) . 
To get the corresponding eigenvalue equation one substi- 
tutes n by n — 1 everywhere in Eq.(25). Everything pro- 
ceeds like in section 0, including the fixed point equation 
E(u>) = uj. But one has achieved much more: Eq.(p7|) is a 
recursion relation which holds for all 1 < n < N\ 

The method of iterated resolvents P7|,p8|j29(] is applica- 



ble to any many-body theory. But due to the Fock-space 
structure as displayed in Table ^ it is particularly easy 
to apply it to gauge theory. Let us demonstrate that in a 
stepwise procedure. For K = 1 the Fock space has only 
one Fock-space sector, and the effective Hamiltonian is 
simply the kinetic energy, i.e. Hi = T\. For K = 2 one 
finds in Table |l| that the Fock space has Nk = 4 sectors. 
By definition, the last sector contains only the diagonal 
kinetic energy, thus H4 — T4. Its resolvent G^ui) is cal- 
culated trivially. Then P3 can be constructed, and thus 
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(73(0;) by a matrix inversion, followed by H 2 and finally 
Hi. Dropping lo in the notation of the propagators, for 
simplicity, one gets consecutively 



Hi = T 4 , (28) 

H 3 = T 3 + VdV , (29) 

H 2 = T 2 + VG 3 V , (30) 

#1 = Tt + VG 3 V + VG 3 VG 2 VG 3 V. (31) 



The effective Hamiltonian Hi with its two terms is re- 
markably simple particularly when compared with the in- 
finite series of the Tamm-Dancoff approach in Eq.(Q). 
The method of iterated resolvents makes efficient use of 
the sparseness of the gauge field Hamiltonian and its many 
zero block matrices. 

To get the effective Hamiltonian(s) for harmonic res- 
olution^) K — 3,4,... is not repeated here explicitly. 
Important is the general feature that the effective sector 
Hamiltonians arc separable in the kinetic energies T and 
the effective interactions U(u)), 

H n {uo)=T n + U n {u). (32) 

Important is also that the effective Hamiltonians in the 
lower sectors become independent of K. The transition to 
the continuum limit K — > 00 is then rather trivial and will 
hence forward be assumed. For the sake of future applica- 
tions the effective interaction was calculated for the first 
12 sectors. Grouping them in a different order, one finds 
with the short-hand notation of Table |l| 

Ui = VG 3 V + VG 3 VG 2 VG 3 V, (33) 

U 3 = VG 6 V + VG 6 VG 5 VG 6 V + VG 4 V, (34) 

U 6 = VG W V + VG W VG 9 VG W V + VG 7 V, (35) 

Uio = VGi 5 V + VGi 5 VG u VGi 5 V + VG U V, (36) 

for the sectors with one gg-pair, i.e. for qq, qqg, qqgg, 
and qqggg, respectively. The effective interactions in the 
sectors with two qq-pairs are: 

Ui = VG 7 V + VG 7 VG 6 VG 7 V, (37) 
U 7 = VGnV + VGnVGioVGnV + VG 8 V, (38) 
U u = VG W V + VGi 6 VGi 5 VG w V + VG 12 V, (39) 

for qqqq, qqqqg, and qqqqgg. Those with three gg-pairs 
have 

U a = VG12V + VG 12 VGuVG 12 V, (40) 
U12 = VGnV + VGi 7 VG m VGi 2 V + VGi 3 V, (41) 

for qqqqqq and qqqqqqg. The structure of the effective in- 
teraction in the pure glue sectors is different: 

U gg = U 2 = VG 3 V + VG 5 V, (42) 
U gg g = U 5 = VG 6 V + VG 9 V , (43) 
Ugg gg = U 9 = VG W V + VGnV . (44) 

One notes finally that Eqs. d28|)-(|3l]) can be reconstructed 
easily from the above relations by setting formally to zero 
all propagators G n with n > 5. 




1 3 2 3 1 




1 3 11 3 



Ui,i u 1>2 



Fig. 1. The three graphs of the effective interaction in the 
gg-space. — The lower two graphs correspond to the chain 
U = VG 3 V, the upper corresponds to U tt = VG3VG2VG3V. 
Propagator boxes are represented by vertical dashed lines, with 
the subscript 'n' refering to the sector numbers. 

Due to the peculiar structure of the gauge field Hamil- 
tonian in Table [j], the vertex interaction appears in the 
effective interactions only in even pairs, typically in the 
combination VGV. It is then plausible that by simply sub- 
stituting 

VGV + VGV (45) 

one restores the instantaneous interaction W which was 
omitted thus far. This rule was checked explicitly in rather 
laborious calculations in (29) . It finds its counterpart in the 
rules for perturbative diagrams []lo[ , where every intrinsic 
'dynamic' line must be supplemented with the correspond- 
ing 'instantaneous' line order by order in perturbation the- 
ory. 

The most important result of this section is that gauge 
theory particularly QCD has only two structurally differ- 
ent contributions to the effective interaction in the qq- 
space, see Eq.(|33|). The effective one-gluon exchange 

U = VG 3 V (46) 

conserves the flavor along the quark line and describes all 
fine and hyperfine interactions. As illustrated in Fig. [l]thc 
vertex interaction V creates a gluon and scatters the sys- 
tem virtually into the qq 5-space. As indicated in the figure 
by the vertical line with subscript '3', the three particles 
propagate there under impact of the full Hamiltonian be- 
fore the gluon is absorbed. The gluon can be absorbed 
either by the antiquark or by the quark. If it is absorbed 
by the quark, it contributes to the effective quark mass 
777. The second term in Eq.(|3^), the effective two-gluon- 
annihilation interaction, 

U a = VG 3 VG 2 VG 3 V, (47) 

is represented by the graph {J 10 in Fig. |l|. The virtual 
annihilation of the qg-pair into two gluons can generate 
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an interaction between different quark ffavors. As a net 
result the effective interaction scatters a quark with hclic- 
ity \ q and four- momentum p = (xP + , xP± + ki,p~) into 
a state with X' q and four-momentum p' = (x' P + , x'P± + 

k^,p' ), and correspondingly the antiquark. In the con- 
tinuum limit, the resolvents are replaced by propagators 
and the eigenvalue problem H e g\ip) = M 2 \ip) becomes 
again an integral equation, but a very simple one in only 
three continuous variables (x, kj_). Its structure is rather 
transparent: 



Af?(aj,kj.;A gs A ff |Vi> 



(x,k±; X q ,Xq\ipi) 
fdx'd 2 ^ 0(x',k' x ){x,k x ;\ q ,\q\VG 3 V + 

VG 3 VG 2 VG 3 V\x\k^X' q ,X' q ) (z',ki; X' q , A^i). (48) 

The domain of integration is set by the sharp cut-off func- 
tion given in Eq.jT^). The eigenvalues refer to the in- 
variant mass Mi of a physical state. The wavefunction 
(x, kj_; A g , Xg\^i) gives the probability amplitudes for find- 
ing in the qq-state a flavored quark with momentum frac- 
tion x, intrinsic transverse momentum kj_ and helicity X q , 
and correspondingly an anti-quark with l—x, —k± and Ag. 
Both the mass and the wave- functions are boost-invariant. 



5 Perturbation theory in medium 

Here seems to be a problem: For calculating G 3 one needs 
Gq, G5 and G4, see Eq.( 



G10, G 9 and G7, see Eqj 



), for calculating Gq one needs 
(34 ), and so on. This property cor- 
responds to some extent the infinite series of the Tamm- 
Dancoff approach, see Eq.(Q). But the method of iter- 
ated resolvents resumes them systematically without dou- 
ble counting. In the next section will be shown how the 
hierarchy can be broken in a rather effective way. That 
final step will be comparatively simple if one has analyzed 
the propagators for the sectors with one qq-pair and arbi- 
trarily many gluons, as follows next. 

Consider first the case with one gluon as given by 
Eq. ( |34] ) . The corresponding diagrams can be grouped into 
two topologically distinct classes, displayed in Figs. Ji] and 
||. Adding one free gluon to the diagrams in Fig. ^ pro- 
duces the diagrams in Fig. |2|. The gluon does not change 
quantum numbers under impact of the interaction and 
acts like a spectator. Therefore, the graphs in Fig. || will 
be referred to as the 'spectator interaction' U 3 . In the 
graphs of Fig. || the gluons are scattered by the interac- 
tion, and correspondingly these graphs will be referred to 
as the 'participant interaction' U 3 . Thus, U 3 = U 3 + U 3 . 

Next consider the propagator G§ in the qq gg-sp&cc. By 
drawing all possible diagrams of Uq, as given in Eq. (p5|) , 
one realizes quickly that most of them topologically belong 
to one of the two classes which one obtains by adding a 
free gluon to the diagrams in Figs. || and ||. The effective 



interaction VGiqV , however, generates also an interaction 
between the two Fock-space gluons in qqgg, by an effec- 
tive one-gluon exchange. Potentially, these diagrams con- 
tribute to a gg-bound-state (a glue-ball), very much like 
the gluon-exchange between the quark and the antiquark 
contributes to a qq-bound-state. By definition, these glue- 
ball diagrams will be included in the spectator interaction 
XJ% . Estimates on their size will be given elsewhere [M . 

The separation into spectators and participants can be 
made in all quark-pair-gluon sector interactions: 



U n = U n + U n , for ri = 3, 6, 10, 15, . 



(49) 



More explicitly, the spectator interactions in the lowest 
sectors with one quark-pair become for example 

U 3 = VGqV + VG 6 VG 5 VG 6 V , (50) 
U 6 = VG 1Q V + VG W VG 9 VG 10 V . (51) 

Correspondingly the participant interactions are 

U 3 = VG4V +VG 6 V , (52) 
U e = VG 7 V +VG W V. (53) 

The same operators appear in both U and U. 

Since the Hamiltonian is additive in spectator and par- 
ticipant interactions, U n can be associated with 



G r . 



1 



cj - T n — U n 
(while G„ = 5—— = = — 



(54) 



The relation to the full resolvent G n = G n + G n U n G n is 
exact. Equivalently, it is written as an infinite series 



G n — G n + G n U n G n + G n U n G n U n G n 



(55) 



The difference to the perturbative Tamm-Dancoff series is 
that the 'unperturbed propagator' Go (a;) in Eq.(p4]) refers 
to the system without interactions while here the 'unper- 
turbed propagator' G n contains the interaction in the well 
defined form of U n . One deals here with a perturbation 
theory in medium. 

The different physics should be emphasized: The sys- 
tem is not scattered into other sectors, it stays in sector 
n. This is reflected in an operator identity, 



(LO-H n )(l-UnGn)=LO-H n 



(56) 



which is obtained straightforwardly from Eqs.(|49|) and 
( B4j ) . The inverse gives R n G n = G„ , with 



Rn 



1 ~ G n U n 



(57) 



With the obvious identity R n G n = G n R n i one ends up 
with G„ = R n G n R n - In all of the above effective inter- 
actions the R n are sandwiched between a quark-pair-glue 
resolvent G and the vertex V, 



VG n V^V R n G n R ] n = V G n V* 



(58) 
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Fig. 2. The three graphs of the spectator interaction in the 
gqg-space. Note the role of the gluon as a spectator. 



Each block matrix V is rectangular and multiplied with a 
square matrix R n according to 

V = VRn, or V ] = R ] n V^ . (59) 

In the sequel we shall suppress again the dagger indicat- 
ing the hermitean conjugate matrix. In Section [?] will be 
shown that R n is essentially diagonal and independent of 
the spin, such that each vertex matrix element is multi- 
plied with a number, actually with a number which de- 
pends on the momentum transfer Q across the vertex. 
Equivalently one replaces the coupling constant g by an 
effective coupling constant g, i.e. 

g^g{Q) = gRn{Q)- (so) 

R n has thus the interpretation of a vertex function. 

One can thus rewrite Eqs.(|35|), ( |50| ) and ( p3l| ) in such a 
form that they are all essentially equal, i.e. 

U e = VG W V + VG w VG g VG 10 V , (61) 
U 3 = VG 6 V + VG 6 VG 5 VG 6 V , (62) 
Ui = VG 3 V + VG 3 VG 2 VG 5 V . (63) 

The effective sector Hamiltonians H n = T n + U n describes 
bound states of one qq-pair with arbitrarily many gluons 
and glue balls. Approximatively, one can relate them to 
each other, see Sec. ||. 

The content of Sects. |] and |B| is exact but rather for- 
mal. To show its usefulness, rigor will be given up in the se- 
quel in favor of transparency and the aim to obtain a sim- 
ple and solvable equation. It should be emphasized here 
already that the content of Sects. |6| and [?] will have to be 
substantiated in future work Bcf|. 



6 The breaking of the propagator hierarchy 

All effective sector Hamiltonians can be diagonalized on 
their own merit. To shape notation, the first few eigenvalue 



1 




I 






(33) 4 


(34) 4 

i 










(35) 6 


(36) g 
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(37) 6 


(38) 6 



Fig. 3. Some six graphs of the participant interaction in the 
gg<?-space. 



equations are written down explicitly. In the gy-space they 
read in Fock space representation 

M%i (?; q\4 1] ) = $>; q\HiW; («'; ?U (1) >- (64) 

For simplicity, the eigenvalues are enumerated by i = 
0, 1, 2, . . . , oo. In the sequel, the lowest eigenvalue is sup- 
posed to satisfy the fix-point equation ui = Mf (u>) and its 
numerical value will be denoted by M 2 . Correspondingly, 
one has in the gg-space 

Ml,, (g;g\^ 2) ) = $>;slW;s') WiS'l^h (65) 
g',§' 

and in the qq g space 

Ml^q-q-g^f) 

= £ (q^gmW^^g^iq'^g'^f). (66) 

i',i',g' 

Knowing the complete sets of eigenfunctions, one can cal- 
culate the exact resolvents in Fock space representation, 
as for example the resolvent in the two gluon sector 

(g;g\G 2 (u)\g';g') 
= Y,^f ) )^^^f ) \9'-r 9 l ). (67) 

j 2;j 

As often in many-body physics, one approximates resol- 
vents by assuming that all eigenvalues are degenerate with 
the lowest state, which here is the glue ball mass M 2 . One 
then applies closure 

(9;g\G2\g';g') = M2 [ M2 J2^9W(^\g';g') (68) 

9 j 
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in order to obtain a diagonal resolvent 



(g;g\G 2 \g';g') = 



1 



M 2 



Ml 



(g;g\g';g') 



(69) 



The propagator in the qqg-sp&ce could be calculated 
by the same procedure, but one can do better. Since the 
gluon is a free particle which moves relative to the qq- 



bound state, the eigenfunction 
product state. The ggg-eigenvalues 



(3)\ - u/,« 



M 



3;j 



M 3 ., s = 



ML 



(i-y) 



W J ) ® \<p. 

7 

y 



is a 



(70) 



can therefore be expressed in terms of the gq-eigenvalues. 
Every qq-bound state is band head for a continuum of 
gluons. With the four-momenta q^ = {yP + , yP± +qj_, Q.g) 
one gets for the lowest bound state 



y z M 2 + qf 



uo - M 3 . = M — M 3 . s - 

y(l -y) 



(71) 



Assuming a degenerate spectrum and performing closure, 
(q;q\g\G 3 \q';q;g') = G 3 (q;q;g)(q;q;g\q';q;g'), (72) 
one gets the delta function (q; q; g\q'; q 1 ; g') multipied with 

G 3 (q;q;g) = Q 2 = (y 2 M 2 + qj 2 )— ^. (73) 

Note that Eqs.(|69|) and ( |73| ) break f/ie hierarchy of the 
iterated resolvents. For calculating the effective interaction 
in the qq-sp&ce only the two resolvents G 3 and Gi are 
needed. Both are written now in closed expressions. The 
whole complication of having resolvents of resolvents is 
replaced by the problem of knowing the eigenvalues M 
and M g ahead of time. Good initial guesses (M ) might 
suffice but can be improved iteratively with M n — > M n+ \. 

The notation in Eq.(f73|) is suggestive for Q 2 being re- 
lated to the single-particle four-momentum transfer along 
the quark line. The free propagator in the qqg-sp&ce can 
be written HI 



G 



1 



qqg, free 



P+(p -pf -q ) (p-p 1 ) 



A2 



(74) 



The single-particle notation refers to Fig. |ij For suffi- 
ciently small y holds (p — p 1 ) 2 — — \y 2 {2m) 2 + qj 2 ]. If one 
substitutes M ~ 2m, which holds to rather good approxi- 
mation, the momentum transfer in Eq.(|74]) is the same as 
in Eq.(|73|). One concludes: In the solution, the interacting 
particles propagate like free particles to a high degree of 
approximation; they just acquire an effective mass. 

Most importantly, instead of having resolvents of re- 
solvents, the hierarchy of iterated resolvents is broken. Ex 
post, this justifies the ui — uj* trick in the Tamm-Dancoff 
work Ull- 

One can restore the exact ggg-propagator by addition 
and subtraction, which gives 

(q;q;g\G3\q';q';g') =G 3 (q;q;g) 

x [(q; q- g\q'; q'; g') - (q; q; g\A\q'; q'; g')] . (75) 



In Eq.(|72|) the operator A was neglected, with 



2 + y{M 2 l -M 2 ) W < sl 



(76) 



Probably, this is the most drastic assumption in this work. 
It suppresses the appearance of the Lamb-shift. 

Note that the above approximations can be controlled 
e posteriori: Once the eigenfunctions in the gg-sector are 
available numerically, one can check whether the exact def- 
inition of the resolvent Eq.(|72]) is peaked like a (5-function, 
with a residue predicted by Eq. (|73|) . 



7 Perturbative analysis of the vortices 

Having established that the propagator G 3 can be approx- 
imated by the free propagator, one can calculate the vertex 
function R 3 straightforwardly. According to Eq.(|57]) one 
has 



R* = 



= 1 + \g 3 U 3 + ^G 3 U 3 G 3 U 3 

1 - G 3 U 3 2 8 



(77) 

It is reasonable to restrict to the first non-trivial term in 
the expansion 

V = VR 3 -V +\ (VG3VG4V + VG 3 VG 6 V) . (78) 

Since the contributions from the vertex functions R4 and 
Re generate terms of higher orders in g, one must set 
i?4 = Rq = 1, in order to be consistent. A diagrammatic 
analysis of Eq.(|7|) un- reveals quickly that it is the famil- 
iar 'vertex correction', which has been calculated in the 
front form by Thorn |37| and by Perry [Q. Thorn and 
later Perry, however, were motivated primarily by asymp- 
totic freedom 'on the light-cone' and thus have emphasized 
the behaviour at large momentum transfer Q — > 00. In 
the present context one is interested in the opposite limit 
Q — > 0. In a bound state equation such as Eq.(^) the 
region near the Coulomb singularity Q ~ is very impor- 
tant. Raufeisen |3^] has recalculated therefore all about 
22 diagrams of the vertex correction in light-cone pertur- 
bation theory, see for example jllj . The salient features of 
these calculations will be exposed here at hand of the two 
vacuum polarization diagrams in Figs. I and |. 

Next let us calculate the free propagator in the qq qq- 
sector, i.e. 



Ga 



1 



(p -p' - q q - q q )P+ 



(79) 



Parameterizing in Fig. ^ the intermediate particle mo- 
menta as 

q% = (yzP+, yzP ± + zqi + q~) , q% = 
(y(l - z)P+, y(l - z)P ± + (1 - z)qx - l ± ,qr) , (80) 
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Fig. 4. The qq vacuum polarization graph. Fig. 5. The gg vacuum polarization graph. 



and using the four-momentum transfer along the quark 
line Q as defined in Eq.(|73]) the propagator becomes 



z(l-z) 



P+ z(l - z)Q 2 + l 2 +mj 
Evaluating the spinor summations gives for rif flavors 



(81) 



/=1 



VG 3 VG 4 V 

1 - 2z(l - z 

G(z,l x ) 



dz dli 



ti+mj + z(l-z)Q 2 



(82) 



Note the spin-independent factor to the matrix element V. 
The cut-off function was defined in Eq. (|l2|) . The term 
in the square brackets are cancelled by other diagrams 
and will be disregarded in the sequel. Integrating over l±_ 
yields straightforwardly 

G 3 VG 4 V = - V-^ 



dz ( : i 



In- 



A 2 + Ami 



(83) 



z(l-z) 



The integral over z is non-trivial. Since the logarithm is a 
very weak function of its arguments, the term m 2 /z(l — z) 



is replaced by its maximum at z = h, thus 



g 3 vg 4 v =-y^ in (i 



A 2 



Am 2 



Q 2 



(84) 



The corresponding steps for the propagator in the qqgg- 
sector (Fig. |) give with G 6 = -q~ -q~)P + ] 



r - H 9 

<~rp, — 7T~ 



z(l-z) 



P+ z(l - z)Q 2 + lj 



(85) 



For the purpose of regularization the gluon is endorsed 
with a small regulator mass m g , see Sec. ^. The spinor 



sums yield 



VGsVGqV = V- 



2n 



dz dl i 



z(l 



0(z,li_) 



l± + m 2 g + z(l - z)Q 2 



(86) 



Dropping the term in the square bracket, and performing 
the same approximation as above gives 



G 3 VG 6 V 



llan c 
12tt 



In 1 



A 2 



Ami + Q 2 



(87) 



As a net result, the coupling constant g has to be replaced 
by the effective coupling constant g, i.e. g — ► g with 



g = gR 3 = g[l 



1 



HQ) 



see also | p9| . The function b(Q) depends on the momentum 
transfer (and on the cutoff), i.e. 



HQ) 



lln e 
12tt 



In 1 + 



A 2 



Ami + Q 2 



— Vln (l 

6tt ^ V 



.4 2 



Am 2 



Q 2 



(89) 



and generates therefore a new interaction. Finally, the ef- 
fective quark masses become according to the diagram 
U 1<2 in Fig. 



—2 2 2 " "c - 

m f = m f + m f — 

1 1 1 7r 2n, 



1, /L 2 

— In — - 

mi 



A similar diagram for the effective gluon mass gives 

2 2 a 2 i ^ 2 

m„ = m„ > m r In 1 J 



4m 2 



(90) 



(91) 



Both are obtained by light-cone perturbation theory [[IT] 
combined with vertex regularization. 

The above considerations are perturbative, as men- 
tioned. What happens if one substitutes the free propa- 
gators in Eqs.(|8l|) and ( |S5| ) by the non-perturbative prop- 
agators G4 = {UJ - Hqqqqy 1 and G 6 = (W - Hqqgg)- 1 , 
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at least in an approximate fashion? — There are addi- 
tional graphs. In the fermion loop of vacuum polarization 
appear two graphs in addition to Fig. ^. In one of them, 
a gluon is emitted and absorbed on the same quark line 
which changes the bare quark mass rrif into the physical 
quark mass m/. In the other graph, the gluon is emit- 
ted from the quark and absorbed by the anti-quark which 
represents an interaction. In consequence one has a bound 
state with a physical mass scale fif. We replace therefore 



2m 



2m/ 



(92) 



Similar considerations hold for the gluon loop in Fig. ^| 
and lead to the substitution 



2m„ 



2m„ 



AV 



(93) 



Both n g and /if are interpreted as physical mass scales. 
The physical gluon mass m g vanishes of course due to 
gauge invariance. This is not in conflict with f.e. Corn- 
wall's suggestion of a finite effective gluon mass since 
one can define (m q ) c s = j-J.g/2. Estimates are given in 



Eq.( |l03| ), and discussed in Sec. ^^n conclusion, for suffi- 
ciently large A, one replaces Eq.(8 




lln / A 2 



6?r ^ 



A 2 



M/ + Q 2 



For Q ^> (fi g , [if) this becomes 
•yl 2 



(94) 



6(Q) = boln, 



with b^ 3 !^!, (95) 
LZir 



for n c = 3. It coincides with the perturbative expressions 
of Thorn and of Perry. The notation is inspired by the 
relation to the familiar beta- function [Ell . 



8 Renormalization 

Now, all the pieces are together which are needed for a 
further discussion of the effective interaction as defined 
in Eq.(|4"8|). In the present first assault the flavor chang- 
ing interaction U a is disregarded. In the remainder U, 



the instantaneous interaction is restored by U„ 



V G q qgV, with G q qg given by Eq.([73j). Both terms contain 
a non-integrable singularity which however cancel each 
other such that only the integrable Coulomb singularity 
Q~ 2 remains, see for example jll). When substituting in 
both of them g — > p one gets from Eq.(f48|) the final one- 
body equation: 



M 2 (a;,k ± ;A„,A, 



m 2 



m| + k 2 
l-x 



,kuAo,A 5 |^) 



-3^2 E dx'd 2 ^G(x>,k' ± ) 
g(g) 

Q 2 



(X q ,X s \S(Q)\X' q ,\L) (x',^;X' q ,X' q \^). (96) 



The effective coupling constant a(Q) = 7? 2 (Q)/47r has thus 
far been give n in Eq. ( |88| ) ; its renormalized version will be 
given in Eq.( [l08|) , below. The cut-off function sets the 
domain of integration and was defined in Eq.([l2|). The 
spinor factor S(Q) represents the familiar current-current 
coupling which describes all fine and hypcrfine interactions 



(Xq,Xq\S(Q)\X , q ,X / q) = 

\u{k q ,Xq)^u{k' q ,X'q)\ [u(k q , Xg)-f^u{k' qi X' q )] 
VxX 7 - X)(l - X 1 ) 



(97) 



Actually, the precise analysis of the formalism requires to 
use the arithmetic mean from each vertex, i. e. 



a(Q) 

Q 2 



-\l a{Q q )a{Q q ) ( g| + g| 



1 



(98) 



where Q 2 q = -(p q - p' q ) 2 and Q\ = -(p q - p' q ) 2 are the 
momentum transfers along the quark and the antiquark 
line, respectively. Close to the Coulomb singularity, how- 
ever, this does not matter since both are approximately 
equal: Q 2 q ~Q 2 q ~ (m q + rriq) 2 (x - x' f + (k ± - k[) 2 . 

Let us return to Eq. (Q) . The analysis there was incom- 
plete and limited to the lowest perturbative order rather 
than to resume the series to all orders in a as required by 
the definition of the vertex function, Eq.([77|). One can cal- 
culate the second order diagram if one neglects all genuine 
two-loop contributions. More generally, if one neglects all 
n-loop contributions one can calculate the term of n-th 
order by 

(G3VG4V + G 3 VG e V) n ~ (ab(Q)) n , (99) 
and resume the series formally to all orders, i.e. 



g(Q) = 



V 1 - a KQ) 



, thus a(Q) = 



l-ab{Q) 



(100) 



The restriction to the one-loop level is a customary pro- 
cedure and not coupled directly to the size of a. More 
important is the point that the series converges only if 



ab(Q) < 1 , 



(101) 



for all Q a nd a . Let us therefore discuss shortly to which 
extent Eq.(lOl) holds true. The possibility of a 1 is 
disregarded, since considerations should not be limited by 
perturbative considerations. According to Eq.(p4|), b(Q) is 
largest for Q = 0. Defining a weighted average [i by 



bo In 



A' 



33 , (A 2 

In 

12tt 



()7T 



/ 



A- 



-Vln. , 



(102) 



gives fi > A 2 e~^o from Eq.(lOl). This lower limit is 
not necessarily small. An estimate appropriate for char- 
monium (M c ) (with A ~ M c and abo ~ 0.2) rather gives 
fi > 245 MeV. Denote the six experimental masses for the 
neutral scalar mesons with Mf and the estimated glue- 
ball mass with M g ~ 1.4 GeV. The scales [if and [i „ are 
identified with them for simplicity. One gets from Eq. (|102[) 



Ms 



M„ 



M 



thus 



0.91 GeV, (103) 
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independent of A and consistent with the above limit. 

This is as far as one can go with the regulated the- 
ory. Next let us address to the renormalization of the the- 
ory, particularly to the renormalization group equations. 
The eigenvalues of the light-cone Hamiltonian Hlc are 
given by Eq.(p^). Potentially, they depend on A through 
three sourcesi(l) the physical flavor masses to/; (2) the 
(9-function representing the domain of integration; (3) the 
effective coupling constant a(Q) . Only the last one is rel- 
evant, since the flavor masses are renormalization group 
invariants with drrif/dA — 0, and since (at least in a con- 
fining regime) the wavefunctions decay sufficiently fast to 
serve as a natural cut-off. It is thus reasonable to replace 
Eq.(i by 

-^a{a;A)=Q. (104) 

It can be satisfied by a cut-off dependent bare coupling 
constant g = g(A), thus by 



da 
dA 



1 



(l-ab(Q)Y 



da 
dA 



, db 
dA 



= 0. 



(105) 



Since db/dA is completely independent of the mass scales 
fig and /x/ appearing in Eq.(g4J), i.e. 



db(Q) 
dA 



A ' 



(106) 



one gets the differential equation da/dA — — 2a 2 bn/A, as 
usual, including its solution El] 



OLA 



1 



1 + ao&o In 



.i 2 



6 In 



A 2 

TP 

QCD 



(107) 



The renormalization point ao at the scale ylo is arbi- 
trary. The two parameters (ao,ylo) can therefore be com- 
pounded into a single one, into the QCD-scale Aqcd = K 
whi ch is defined by 2qq6q I n Aq/ Aqcd = I- Substituting 
Eq.(|l~07|) back into into Eq.pX) gives 



a(Q) = 

and thus 
a{Q) = 



b In- 



A 2 



HQ) 

br, In 



12tt 



b n In ■ 



A 2 



HQ) 



33 In 



Q< 



2£ln 

/ 



fj>j + Q 



(108) 



All In (A) terms cancel exactly without smallness assump- 
tion. The c ancelation is directly related to the minus sign 
in Eq. (|lOC| ), which in turn is related to the re-summation 
of the series to all orders in a. Using Eq.(pq), the cance- 
lation would hold only perturbatively, for sufficiently small 
a. 

Since a becomes independent of A, one can go to the 
limit A — > oo. In line with the modern interpretation of 



the renormalization group one can integrate Eq.(|9q) un- 
restrictedly over all Q. Note the subtle difference between 
the 'running coupling constant' g^ and the 'effective cou- 
pling constant' ~g{Q). The two are confused often in the 
literature. The Lagrangian does not know about the mo- 
mentum transfer. 

A final word on the complete wavefunction The 
solutions of the integral equation Eq.(96) represent the 
normalized projections of the full eigenfunction onto 
the Fock states \q;q) = b\d\\vac), i.e. tp^ = Un- 
like in other methods this knowledge is sufficient in the 
present method to retrieve all other Fock-space compo- 
nents (n\H/) without solving another eigenvalue problem. 
As an example, let us work that out explicitly, by asking 
for the probability amplitude to find a \gg)- and a \qqg)- 
state in a full eigenstate I'P). 

The key is the upwards recursion relation Eq.(p6| ) . The 
first two equations of the recursive set in Eq.(Eq) are 



(2|!P) = G a <2|fr a |l)<l|!F), 
<3|^) =G 3 (3|ff 3 |l>(l|#>- 



G3(3|i7 3 |2)<2|^) 



(109) 
(110) 



The sector Hamiltonians H n have to be substituted from 
Eqs.(j34|) and (^2|). In taking block matrix elements of 
them, the formal expressions are simplified considerably 
since many of the Hamiltonian blocks in Table [l] are zero. 
One thus gets simply (2|H" 2 |1) = (2\VG 3 V\1) and there- 
fore (2\W) = G 2 VG 3 V(l\W). Substituted into Eq.jl09|) 
this gives (3|<f) = G 3 V{1\&) + G 3 VG 2 VG 3 V(l\W r ). These 
findings are summarized more readable 

{ggW) = (gg\G gg VG qqg V (ill) 
(Q99\¥) = (QQ9\G q qgV\% q ) 
+(qqg\G qqg VG gg VG qqg V (112) 

Correspondingly, one is able to find the Fock-space am- 
plitude in the higher sectors with remarkably little effort. 
Of course, one has to readjust the overall normalization 
of the wavefunction. It should be emphasized that the fi- 
nite number of terms is in strong contrast to the infinite 
perturbative series. Iterated resolvents resume the series 
to all orders in closed form. It also should be emphasized 
that the same approximations as discussed above must be 
made for reasons of consistency, particularly the effective 
coupling constant g(Q) must be used. 



9 Summary and discussion 

The present work is based on Lagrangian gauge field the- 
ory for SU(N) and on its canonical front-form Hamiltonian 
in the light-cone gauge A + = 0. The notational back- 
ground is laid down in Sects. || and ||. All zero modes 
are disregarded, and in consequence the vacuum is trivial. 
Imposing periodic boundary conditions the Hamiltonian 
is converted to a matrix. The matrix is finite by means 
of Fock-space regularization. All possible divergences are 
controlled by vertex regularization. Being confronted with 
the diagonalization of a finite but large matrix, the prob- 
lem is mapped on a smaller space as was done first by 
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Tamm and by Dancoff in their theory of effective interac- 
tions. Binding effects arise then due to virtual scattering 
into higher Fock-space sectors. The apparent difficulties 
with this approach are related to the infinite perturbative 
series one is forced to work with in practice. A re-analysis 
shows that the main idea can be maintained if one intro- 
duces a hierarchy of effective interactions, in each sector 
of the Hilbert space. The result is an iterative procedure, 
called the method of iterated resolvents. Each resolvent 
is the inverse of an effective sector Hamiltonian which in 
turn is a functional of resolvents in higher sectors. For the 
model case of a finite matrix (DLCQ) the method of iter- 
ated resolvents can be realized and explicitly checked by 
a finite number of successive matrix inversions and multi- 
plications. But even for the continuum it is a well defined 
and exact procedure. The infinite perturbative series of the 
Tamm-Dancoff approach and all their many-body aspects 
are then replaced by a finite number of terms. The effective 
interaction between a quark and an antiquark turns out to 
have only two contributions: The flavor conserving interac- 
tion U and flavor changing interaction U a . Their diagram- 
matic representations look like second order diagrams of 
perturbation theory, but represent a re-summation of per- 
turbative graphs to all orders. Particularly U bears great 
similarity with a perturbative one-gluon-exchange interac- 
tion. As part of the approach the complete wavefunction 
can be reconstructed, by evaluating one component after 
the other in a well defined procedure. Once this is achieved 
one can relax periodic boundary conditions and return to 
the continuum limit. The result is a complete and exact 
theory of the effective interaction between a quark and an 
antiquark. ft can be interpreted as the genuine interac- 
tion in a constituent quark model. The body of this work 
is found in Sects. || and |[ 

In Sects. ^ and essentially four assumptions are made 
for the sake of transparency which shall be summarized in 
brief. — The neglect of the operator A in Eq. ((75]) is proba- 
bly the most drastic assumption. It prevents a straightfor- 
ward calculation of the Lamb shift and reminds to propo- 
sition 13 in Feynman's famous article (43[ |. It is here where 
statistical considerations |Q can perhaps be implemented 
in the future. The next severe approximation resides in 
Eq. ([Tp) , where the eigenvalue M appears in the propaga- 
tor. If one substitutes M ~ 2m, the propagator in the qqg 
space behaves like the propagator of a free particle with an 
effective mass m (rather than with the bare mass m), to a 
high degree of approximation. This step is important since 
it 'breaks the hierarchy': The propagator G q q g need not 
be obtained from an iteration procedure. In a similar way 
one gets a simple approximation for G gg . Only these two 
propagators occur explicitly in the expression for the effec- 
tive interaction. By the nature of the approximation both 
of them become independent of the starting point energy 
uj. The third and fourth assumption resides in the vertex 
function, which contains all many-body effects to arbitrary 
order of the cou pling constant. The vertex function was 
evaluated in Eq. ( 100| ) up to the one- loop level, as done of- 
ten in applications of gauge theory. - The following point 
should be emphasized: Simplifying assumptions are made 



here only after having found the general structure. They 
are therefore unlikely to violate fundamental symmetries 
like gauge and Lorcntz-invariance. Usually one proceeds in 
the reverse order |i^ , |l^ , p7|Jl^ , [l^ , p6| , ^llj3 2j , |34| , [35| : One first 
truncates and then develops the formalism. 

Taking all together in Sect. || one arrives at a compar- 
atively simple integral equation in the variables of a sin- 
gle quark, the one-body equation (|9q) . Its kernel contains 



the effective coupling constant a. It is defined in Eq. (|108| ) 
and accumulates approximatively all many-body aspects. 
It depends on the four-momentum transfer Q along the 
quark line and on the QCD scale Aqcd = k. The latter 
arises by renormalization and must be determined by ex- 
periment. Beyond that, ~a(Q) has a rich parametric struc- 
ture depending on one gluonic and six fermionic mass 
scales, \i a and fij, respectively. They are in principle cal- 
culable Q, but as explained in the context of Eq. (f73|) , 
the present formalism needs an initial input guess which 
later can be improved iteratively and self-consistently. At 
present \i g and \j,f a re ta ken as external parameters. The 
values quoted in Eq. ( 103 ) look like reasonable first guesses. 



But one can view [i g and /i/ also as external param- 
eters, which are fixed subject to describe experiments or 
other theories. This opens a broad avenue for QCD-based 
phenomcnological applications. One knows where the as- 
sumptions have been made in the present formalism, and 
can relax them in subsequent work. 

For example, one can choose \x g = jif = K to get 
&(Q) = fro hi (1 + Q 2 /k 2 ). The numerical value of k = 
0.19 GeV gives a reasonable fit to the empirical masses 
of the scalar and vector mesons [P0[ . This form was intro- 
duced by long ago Richardson |45|| to interpolate between 
asymptotic freedom and infrared slavery, see also (42). It 
generates a Coulomb plus a linear potential in configu- 
ration space, and thus produces confinement. "There is 
a pleasure in recognizing old things from a new point of 
view" H) . Unfortunately the Richardson coupling has the 
unphysical aspect of tending to infinity for vanishi ng fo ur- 
momentum transfer. But the full expression Eq.( |10q ) is 
finite for Q = since S(0) _1 = &ohi(^ 2 /K 2 ), with a nu- 
merical value /i = 0.91 GeV as given in Eq.( |l03| ). The 
parametrization 



l 



6nln^ 



(113) 



is simple and suggestive as an approximation. It has been 
used repeatedly in the past, see f.e. p2| , as reviewed in 
p6| . Brodsky et al |Q have fixed the parameters by a 
fit to the non-relativistic heavy-quark lattice data || and 
get k — 0.16 GeV and /i = 0.87 GeV. The (perhaps ac- 
cidental) coincidence with the above quoted numbers is 
amazing. 

It would be more than interesting to solve the inte- 
gral equation (|96|), or approximations thereof, with various 
parametrization of a(Q), or at least calculate the poten- 
tial in configuration space. The work is under way but for 
obvious reasons must be dealt with separately |4(J. The 
investigation of chiral symmetry (breaking) must also be 
postponed for the future. In the present analysis it can- 
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not be discussed since the flavor changing interaction U a 
was disregarded waiting to be calculated. It also should 
be interesting to apply the method of iterated resolvents 
to conventional many-body problems, and to re-analyze 
QED. 

To summarize the present work in short one can state 
that the effective potential between the constituent quarks 
in a meson has been derived from the bare Lagrangian for 
(actually any) SU(N) gauge held theory. It is the hrst time 
that this was possible within a self-consistent treatment. 
The approach is based a novel technology and has some 
virtues, among them: 

• The minimum number of physical degrees of freedom 
are used because of the light-cone gauge; 

• All Lagrangian symmetries are preserved; 

• Rcnormalization is carried out explicitly; 

• Fermions are treated exactly; thus no fermion doubling; 

• The final one-body equation is solvable. 
A lot of work remains to be done. 
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